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ABSTRACT 

We present the results of wide integral-field near-infrared (1.0-1.8 /um) spectroscopic observations 
of the southeastern shell of the young core-collapse supernova remnant (SNR) Gil. 2— 0.3. We first 
construct [Fe II] 1.644 /ini line images of three bright clumps from the obtained spectral image cubes 
and compare them with those of other transitions such as [Fe II] 1.257, [Fe II] 1.534 and He I 1.083 //m 
line images. This allows us to estimate the electron density (~ 4,700-9,400 cm~'^) and extinction 
(Ay ^ 16-20 mag) of the shell, including detailed two-dimensional distribution of the properties in 
the brightest clump, as well as the discovery of a faint high-velocity (~ -440 km s~^) component in 
the clump. Our SNR shock model calculations estimate the preshock number density of ^ 250-500 
cm~^ and shock speed of ^ 80-250 km s~^ in the [Fe II]-emitting region of the SNR. The comparison 
between the observed and modelled radial profiles of the line intensities and their ratios reveals that 
the shell is composed of multiple thin filaments which have been likely formed in episodic mass 
loss processes of a progenitor star. The discovery of the faint high-velocity component supports the 
interpretation that the southeastern shell of Gil. 2— 0.3 is mainly composed of circumstellar material 
with contamination by supernova ejecta and also that its ejected material was expelled primarily in 
the southeast-northwest direction. 

Subject headings: ISM: individual objects (Gil. 2— 0.3) — ISM: supernova remnants — infrared; ISM 
— shock waves 



1. INTRODUCTION 

Near-infrared [Fe II] observation of young core-collapse 
SNRs in the Galactic plane is an invaluable tool for 
understanding many interesting phenomena associated 
with core-collapse supernova explosions. Iron is the end 
product of the stellar nucleosynthetic processes with rela- 
tively small (7.90 eV) ionization potential, which renders 
near-infrared [Fe II] observations particularly useful for 
studying the distribution and kinematics of ejecta from 
deep layers of a progenitor star and those of circumstellar 
material from its later evolutionary phases prior to the 
explosion, along with interstellar shocks, while suffering 
less extinction effect than the visual or X-rays. 

An enlightening illustration is the heavily extinquished 
{Ay ^ 13 mag) young core-collapse SNR Gil. 2— 0.3, 
the likely remains of the historical supernova explosion 
in SN 386 at 5 kpc distance ([Stephenson &: GreenI 
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[200l IBecker et all 119851 : IGreen et all [T98l . of which 
recent [Fe II] observations have revealed se veral perti- 
nent featur es. Accord i ng to the results of I Koo et all 
11007) and lMoon et"aI1 (l2009l) . Gil. 2-0.3 is filled with 
bright near-infrared [Fe II] emission originating from 
both the shocked ejecta and circumstellar material scat- 
tered around the entire remnant. In its inner part, 
knotty [Fe II] emission with significant ( ^ 1,000 km s~^) 
Dopplcr shifts is prominent without any other bright el- 
ements, indicating that the majority of the iron knots 
therein originate deep inside its progenitor. Near its 
northwestern boundary lie knotty filaments of [Fe II] 
emission, reminiscent of bubbly iron ejecta. The most 
conspicuous near- infrared [Fe II] feature in Gil. 2— 0.3, 
however, is its southwestern shell which is in fact the 
brightest near-infrared [Fe II] emission ever detected 
among all Galactic SNRs. The shell is thick and 
pro minent in radio. X-ray and mid-infrared emission 
too llGreen et al.lll988l : iKaspi et al.ll2001l : iPinheiro et al.l 
12011 ). Although all these are suggestive that the south- 
western shell is a rare example of shocked supernova 
ejecta interacting with shocked circumstellar material, 
the origin and nature of this unique feature are uncer- 
tain. 

One efficient way to investigate the [Fe II] emission of 
the southwestern shell of Gil. 2— 0.3 is to obtain broad- 
band two-dimensional integral-field spectra of the entire 
emission which can provide both spectral and spatial in- 
formation simultaneously. In this paper, we report the 
results of such observations using our own instrument 
the Florida Image Sheer for Infrared Cosmology and As- 
trophysics (FISICA). We describe our instrument setup 
and observations in § 2, followed by results of the obser- 
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vations (§ 3) and discussions focused on shock modelling 
of multiple near- infrared [Fe II] line emission (§4). We 
summarize our results in § 5. 

2. INSTRUMENT SETUP AND OBSERVATIONS 

The instrument setup used for the observations we re- 
port in this paper is a combination of Florida Multi- 
object Imaging Near-infrared Grism Observational Spec- 
trometer (FLAMINGOS) and FISICA. The former is a 
facility instrument at the Kitt Peak National Observa- 
tory Mayall 4 m telescope. One unique feature of this 
instrument is a small multi-object spectrograph (MOS) 
dewar, separated from its m ain dewar, ded icated for cy- 
cling multi-object slit masks (|Elston.lll998l) . We replaced 
the MOS dewar with FISICA which is an image sheer- 
based integral -field unit composed o f three monolithic 
mirror arrays (|Eikenberrv et all [20061 ) . FISICA was de- 
veloped by some of the authors of this paper as an ac- 
cessory module for FLAMINGOS, and the combination 
of FLAMINGOS and FISICA provides an unprecedent- 
edly large integral- field of view of ^ 16" x 33" in the 
near-infrared regime. This setup is uniquely ideal for 
studying kinematics and chemistry of extended objects 
in the Galactic plane, such as the [Fe II] emission of the 
southwestern shell of Gil. 2— 0.3 as we report here. 

The observations took place on 2006 June 17 and 18 
using the JH-hand grism of FLAMINGOS, which cov- 
ers the ~ 1.0-1.8 fim range populated with many [Fe II] 
emission lines with R ~ 1,300 spectral resolving power. 
We observed three bright [Fe II] clumps (which we call 
Clump 1, 2, and 3) within the southeastern shell of 
Gil. 2 -0.3 originally detected in t he [Fe II] 1-644 /xm 
narrow-band imaging observations ()Koo et al.l 120071 see 
Figure [T]) . We followed each exposure of 300 seconds to- 
ward the clumps with the same exposure toward nearby 
sky without any apparent [Fe II] emission in order to sub- 
tract out sky background emission (including OH lines). 
The total on-source integration times were 1,500, 2,100, 
and 900 seconds for Clump 1, 2, and 3, respectively, 
and the same integration times were used for sky ob- 
servations. After flat fielding, we obtained wavelength 
solutions using HeNeAr arc spectra, and used the stan- 
dard star of SAO 1 87086 (G6V) for telluric calibration 
(jCoelho et al.|[2005f) . The astrometric solutions of the 
obtained data cubes are consistent with the known po- 
sitions of stars located inside the integral fields of the 
clumps. 

3. RESULTS 

Figure [T] presents pure [Fe II] (1.644 /im) line images 
of Clumps 1, 2, and 3, together with those of He I 
(1.083 /im) and [Fe II] (1.257 and 1.534 ^m) of Clump 
1. To construct the line images, we measure the inten- 
sities of the line transitions at every spaxel by fitting 
their profiles with a gaussian function. The [Fe II] (1.644 
Hm) line images of the three clumps show extended (> 
25") elongated structures along the southwestern shell of 
Gil. 2— 0.3. For the brightest Clump 1, in addition to 
the 1.644 /im transition, we also detect significant emis- 
sion of [Fe II] from the transitions at 1.257 and 1.534 /im. 
The images of the latter transitions closely resemble that 
of the former, with correlation coefficients of 0.94 (1.256 
vs. 1.644 fim images) and 0.95 (1.534 vs. 1.644 fim im- 
ages). We also detect He I (1.083 /Ltm) line emission from 



Clump 1. Although faint, its peak position is coincident 
with those of the [Fe II] lines, and it appears to extend in 
the same direction as the extension of the [Fe II] images. 

Figures [2^ and 2b show JH-hand spectra of Clumps 1, 
2, and 3 integrated over the entire region of each clump. 
As in Figure m they are exclusively dominated by [Fe II] 
transitions, and in total we identify 12 [Fe II] transi- 
tions. The only line which is detected other than the 
[Fe II] lines is He I (1.083 /xm) from Clump 1. Table [U 
contains the details of the detected line transitions in- 
cluding their relative intensities. The detected lines do 
not show any apparent Doppler shif t s, wh ich is consis- 
tent with the results of iMoon et al.l (|2009[ ) around this 
area of the SNR, indicating that their motions are mostly 
tangential. However, in the profile of the [Fe II] (1.644 
/im) transition from Clump 1 (Figure [5};), we identify an 
additional velocity component which we call the south- 
eastern blueshifted high- velocity feature (or SE-BHVF). 
The extension of SE-BHVF is ^ 10" and its peak posi- 
tion is located at 15" from the main peak position of 
Clump 1 in the northeastern direction (Figure [ij . Using 
a two-velocity component gaussian fit of the profile, we 
estimate the radial velocity of SE-BHVF to be -440 ± 
100 km s^^, highly blueshifted from the systematic ve - 
locity km s^^) of GU. 2-0.3 (| Green et al.l[T988h . 

The integrated intensity of SE-BHVF is < 2 % of the 
main component. 

Because the two [Fe II] transitions at 1.257 and 
1.644 /im detected in Clump 1 share the same upper 
level, their intrinsic intensity ratio is determined by their 
atomic parameters. We compare the observed ratio with 
the intrinsic ratio o f 1.36 expected from the atomic pa- 
rameters of [Fe II] ()Deb &: HibbertI [2011[ ) and compute 
the extinction toward the source using the ex tinction 
curve of grains in the Milky Way with Ry = 3.1 (jPraind 
I2003f ). Figure [3)3 shows a visual extinction (Ay) map of 
Clump 1 which presents two-dimensional distribution of 
the extinction in the clump. Its mean extinction is Ay = 
16 ± 1 mag, but varies somewhat significantly within the 
clump between 12 and 20 mag and increases in the south- 
western directions where Clumps 2 and 3 are located. For 
the extinctions toward Clumps 2 and 3, we only calculate 
their average extinctions of Ay = 18 ± 1 (Clump 2) and 
20 ± 1 (Clump 3) mag due to relatively weak intensities 
of their [Fe II] 1.257 /im transitions. The internal ex- 
tinction distribution within Clump 1 and the increased 
extinctions toward Clumps 2 and 3 show a tendency that 
the overall extinction toward Gil. 2— 0.3 in this area in- 
creases along the southern direction, indicating potential 
existence of a foreground or circumstellar gas cloud there. 
The extinction toward Clump 1 is consistent with a SN 
at 5 kpc being visible to the naked eye in AD 386. The 
average extinction (Ay — 16 mag) of Clump 1 is higher 
than the extinction (Ay — 13 mag) pr eviously estimated 
at the intensity peak position (|Koo et al. 2007), although 
their observed transition ratios of [Fe II] 1.257 to 1.644 
/im are similar within errors (0.303 ± 0.019 for the former 
and 0.314 ± 0.010 for the latter), simply due to the differ- 
ent atomic parameters and resulti ng intrinsic intensities 
used to estimate t he extinctions (jPeb fc HibbertI [20111 : 
IQuinet et anil996D . 

We also use the distribution of the transition ratio of 
[Fe II] 1.533 to 1.644 /tm to obtain two-dimensional dis- 
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tribution of the electron density in Clump 1 (Figure [3J:). 
We solve the rate equations of 16 levels of Fe"*" using the 
transitio n probabilities and colli sion strengths in the lit- 
erature ()Deb fc Hibbert 2011; Ramsbottom et al.l [20071 ) 
at the temperature of 5,000 K. As in Figure [5J;, the dis- 
tribution of the electron density in Clump 1 is totally 
different from that of the [Fe II] 1.644 line intensities 
with their peaks located at different positions. The aver- 
age electron density of Clump 1 is 9,400 ± 2,100 cm~^, 
while those of the other two clumps are 8,900 ± 2,800 and 
4,700 ± 1,000 cm^'^ for Clumps 2 and 3, respectively. 

4. DISCUSSION 

The near-infrared JiJ-band spectra of the observed 
clumps are dominated by [Fe II] lines - only weak emis- 
sion of He I is detected other than [Fe II] lines in Figure[21 
One notable feature is the lack of correlation between the 
distribution of the [Fe II] line intensities and that of the 
electron number density (Figure [3]), which indicates that 
the [Fe II] 1.534 to 1.644 ^m ratio distribution is largely 
independent of the line intensity distribution. In order to 
investigate this more thoroughly, we compare the radial 
distribution of the two quantities, i.e., [Fe II] 1.644 fim 
line intensity and the [Fe II] 1.534 to 1.644 /im line ratio, 
along the peak position of Clump 1 (Figure ^ by apply- 
ing boxcar average of 5". Figure 2] shows that the radial 
distributions of the two quantities are clearly dissimilar: 
the [Fe II] 1.644 /im line intensity peaks at the middle of 
the clump and decays both inward and outward with ^ 
4" FWHM, while the [Fe II] 1.534 to 1.644 ^m ratio dis- 
tribution is almost flat with a tendency of slight increase 
toward the outer edge of the clump. 

We model the radial profiles of the [Fe II] 1.534 and 
1.644 /im line intensities and their ratio for a SNR shock 
penetrating into dense ambient medium using a shock 
co de (Figure HI). T he code that we u se was developed 
bv iRavmondl ()1979D and improved bv I Cox fc Ravmondl 
(|1985f ). We have updated the atomic par ameters of [Fe II] 
forbidden lines, using recent va l ues fro m lDeb fc Hibberg 
((2011J) and Rams bot t om et al.l (|2007f ) for the radiative 
transition rates and the coUisional strengths, respec- 
tively. We calculate the level populations of the lowest 
16 levels in four terms (a^D, a'^F, a^D, and a^P) of 
Fe"*" and the intensities of the lines resulting from the 
transitions among them. The line emissivities of a fluid 
element are calculated following its trajectory from the 
shock front until it cools down to 1,000 K, and are used to 
compute the line intensities and their ratios. We assume 
that incoming H atoms are fully ionized by radiative pre- 
cursor. The code, however, does not calculate the emis- 
sion from radiative precusor itself. This is acceptable for 
[Fc II] lines because Fe atoms in the pre-shock gas may be 
mostly locked up to grains, which are liberated when the 
grains are destroyed in the post-shock flow. The [Fe II] 
1.534 and 1.644 nm line brightnesses and their ratio are 
functions of shock speed (vs), preshock density (no); and 
preshock magnetic field strength {Bq) perpendicular to 
the shock front in general. We will present some gen- 
eral characteristics of the [Fe II] emission from radiative 
shocks using the shock code in the forthcoming paper 
(Koo et al. in preparation). Overall the preshock den- 
sity is well constrained by the [Fe II] 1.534 to 1.644 /xm 
ratio. 

The average [Fe II] 1.534 to 1.644 ratio of Clump 



1 is 0.192 ±0.012 (Table[T]), indicating the average post- 
shock electron density of ^ 1 x 10^ cm~'^ (see § 3). 
This value should be much higher than the uncompressed 
preshock density for a radiative shock. For shock speeds 
of 80-250 km s'-^ and Bq = 10 /iG, the observed ratio 
implies no = 250-500 cm~'^. For _Bo = 1 /iC, it becomes 
no = 150-300 cm~^ for the same shock speed range. We 
adopt 10 /iG, which is the median value o f the maximum 
magn etic field strength for diffuse clouds (jCrutcher et al.l 
|2010[) . We take the shock speed to explain the observed 
line emission: if the velocity is considerably lower, e.g., 
50 km s~^, we may expect to see H2 emission from 
molecular shocks instead of bright [Fe II] emission from 
atomic radiative shocks. If the velocity is siginificantly 
larger, e.g., > 250 km s~^, the expected brightness nor- 
mal to the shock front becomes brighter than the ob- 
served one, i.e., 1-2 xlO~^ ergs cm~^ s~^ sr~^, which is 
th e intensity of the b right structure of Clump 1 obtained 
bv lKoo et al.l (|2007D after being smoothed to the beam- 
size (2") of the current observation and dereddended 
using Ay — 16 mag. Note that the observed bright- 
ness is the integrated intensity along the line of sight 
and also beam-diluted. The brightness along the line of 
sight is much larger than the brightness normal to the 
shock front in general, so that it is not straightforward 
to give more stringent constraints without more compli- 
cated and detailed shock model calculations. 

We perform a shock model calculation of the observed 
radial profiles assuming solar abundance ([Asplund et al.l 
l2009f ). In Figure m the red-solid and red-dashed lines 
represent the normalized [Fe II] 1.644 iim intensity and 
the [Fe II] 1.534 to 1.644 /im line ratio, respectively, for a 
spherically symmetric shock of 130 km propagating 
into an ambient medium of no = 300 cm^'^. The peak 
intensity of [Fe II] 1.644 /im line profile is 8 x 10^"^ ergs 
cm^^ s^^ sr^^. Note that these are profiles integrated 
along the line of sight and beam-diluted. For compari- 
son, we also show the normalized [Fe II] 1.644 /im inten- 
sity profile (black-dotted line) which is free of the beam 
dilution effect, where we can identify a sharp intensity 
contrast caused by a thin shock front. As in Figure |31 
our shock model fits the [Fe II] 1.534 to 1.644 /im ratio 
very well, whereas it does the [Fe II] 1.644 fim intensity 
profile somewhat poorly. The observed intensity (filled 
circles) is much fainter than the modelled one (red-solid 
line) at i? < 0.95 and the width of the observed shell 
is much greater than what is modelled. These discrep- 
ancies imply that the structure is not spherically sym- 
metric, with local enhancement of high-density medium 
in the region, and also that there may be multiple fila- 
ments inside the beam. We can fit the observed intensity 
profile and the peak brightness readily using multiple fil- 
aments of finite sizes. For example, the green-solid line in 
Figure |4] shows the modelled intensity profile when there 
are three filaments of 1.7 pc size separated by 0.03-0.04 
pc from each other, and it is in a good agreement with 
the observed values. The current model is probably not 
a unique model for the radial profile, but multiple fila- 
ments seem to be unavoidable. A possible scenario for 
the origin of such multiple filaments is an episodic mass 
loss from the progenitor of Gil. 2— 0.3 to the southeast- 
northwest direction. The similar distribution of He I in 
Clump 1 (Figure [T]) also agrees with the circumstellar 
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origin of the main [Fe II] emission of the clump, exclud- 
ing SE-BHVF (see below). The dereddended line ratio 
of He I 1.083 to [Fe II] 1.644 /um is ~ 0.50, which is 
comparable to ~ 0.59 previously observed in RCW 103 
(lOliva et al.l[T990l) . 

The three [Fe II] clumps that we observe here are sur- 
rounded by hot X-ray gas (Figured]). The density of the 
X-ray emitting gas is 13 cm~^, if we simply divide its 
ionization time scale rip t = 6.7 x 10^^ s cm ~'^ by the 
SNR age of 1,600 years ([Roberts et al.l 12001 . The tem- 
perature of the X-ray emitting gas is 0.58 keV, so that its 
pressure becomes p/ks — 1.8 x 10^ cm~^ K, where /cs 
is the Boltzmann constant. For comparison, the average 
electron densities of clumps are 4,700-9,400 cm"'^. Us- 
ing the characteristics temperature of 5,000 K, we get the 
ionized gas pressure (rie + np) T = 0.47-0.94 x 10^ cm~'^ 
K, which is about 1/4 to 1/2 of the pressure from the 
X-rays, indicating that the ionization fra ction of [Fe II] 
emitting gas is low fe.g. JOliva et al.|[i989() . The total gas 
pressure, neglecting He, is p/ks = (1 + fe)/ (2/e) x rieTe, 
where fe, is the mean ionization fraction. This gives 
fe = 0.14-0.34, consistent with our previous results from 
the shock model calculations 

The southeastern shell of Gil. 2— 0.3 has been detected 
by Spitzer mid- infrared IRAQ 3-8 /xm and MIPS 24 
jxm i mages (|Leel 120051: [Reach et al.|[200l IPinheiro et al.l 
120111) . Comparing the distributions of the mid- infrared 
images with those of the ne ar-infrared [Fe II] and H2 
line images (|Koo et al.ll2007D . the brightest Clump 1 is 
less pronounced in the mid-infrared images, and a small 
molecular c loud located at y 1' southwest from Clump 1 
(H2-pkl in lKoo et al.ll200"7|) is prominent in the Spitzer 
IRAC images. This is consistent with the classification 
of Gil. 2— 0.3 (H2-pkl) with molecular sh ocks by the ra- 
tios of the IRAC intensities (jReach et al.li2006.) . and the 
Spitzer mid-infrared spectrum, obtained at the position 
of this molecular cloud, shows a wealth of H2 pure ro- 
tational and ionic line s as well as bright dust continuum 
([Andersen et al.l 120111) . Unfortunately, the position of 
the molecular cloud is outside the [Fe II] clumps and not 
covered by our observations. For the Spitzer MIPS 24 
/im emission of Clump 1, although we expect a signifi- 
cant contribution from [Fe II] 24.52 and 25.99 /xm lines, 
it is not straightforward to determine the relative contri- 
butions from these ionic lines and dust continuum, due 
to various uncertainties in comparing the model spec- 
trum to the observed brightness. Future imaging spec- 
troscopic observations are needed for detailed study of 
the mid-infrared emission. 

An interesting finding of our observations is the detec- 
tion of SE-BHVF which is Doppler shifted by -440 ± 
100 km s^^ around the northeastern corner of Clump 1 
(Figures [1] and H:). Assuming that SE-BHVF is lo- 
cated at the boundary of 2' radius of the SNR shell 
with its projected distance of 1.7' from the center, its de- 
projected expanding velocity is ^ —900 km s~^ with re- 
spect to the -|-45 km s~^ systematic velocity of Gil. 2— 0.3 
(jGreen et al.l 119881 ). This expanding velocity is compa- 
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rable to the velocities of the highly blueshifted —1000 
km s~^) ejecta kn ots previously dete cted around the cen- 
ter of Gil. 2-0. 3 (jMoon et al.H2009t ). which is suggestive 
that SE-BHVF has the same origin as the central knots. 
If Gil. 2— 0.3 is inde ed the remnant of the supe rnova ex- 
plosion in AD 386 (jClark fc StephensonI [19771) . the ob- 
served ejecta velocities indicate that they have been con- 
siderably decelerated. The existence of SE-BHVF in 
Clump 1 of the southeastern shell is consistent with the 
previous interpretations that the SN ejecta in Gil. 2— 0.3 
was expelled mainly in the southeast-northwest direc- 
tion, given the existence of the northwester n filaments 
which are most likely of the ejecta origin (|Koo et al.l 
l20(T7l: IMoon eTallllliM . 



5. CONCLUSION 

We perform near-infrared (1.0-1.8 /xm) integral-field 
spectroscopic observations of the southeastern shell 
of the young core-collapse SNR Gil. 2— 0.3 using our 
own instrument setup that provides an unprecedented 
integral-field size in the wavebands. The shell is known to 
have brightest near-infrared [Fe II] emission of all Galac- 
tic SNRs, and is presumed to be a unique mixture of 
progenitor's circumstellar material and supernova ejecta. 
Our data cube populated with multiple [Fe II] transitions 
makes it possible to study the detailed two-dimensional 
distributions of the line intensities, velocities, electron 
densities, and extinction. Based on this information, we 
compare the observed results with those predicted by 
SNR shock models and obtain shock parameters in this 
region of the SNR. The observed properties of the bright- 
est clump can be explained by multiple shocked thin fila- 
ments of circumstellar material. This, together with the 
discovery of the faint high-velocity component, suggests 
that the southeastern shell of Gil. 2— 0.3 is indeed a mix- 
ture of circumstellar material and ejecta, with the former 
being the dominant component. Future studies with im- 
proved sensitivity and resolutions (both spatial and spec- 
tral) are needed to understand how the two components 
are interacting. 
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Fig. 1. — Left: Near-infrared line images of three bright [Fe II] clumps within the southeastern shell of Gil. 2— 0.3. Prom the top, the 
gray-scale images are He I 1.083, [Fe II] 1.257, [Fe II] 1.534, and [Fe II] 1.644 fiicn line images of Clump 1, and [Fe II] 1.644 fim line images 
of Clump 2 and Clump 3. Gray-scale units are normalized to have peak intensities of 100 for the [Fc II] 1.644 /im line image of Clump 1. 
Contours on the [Fe II] 1.644 fim line image of Clump 1 indi cate the distribu tion of the SE-BHVF. Right: [Fe II] 1.644 /im narrow-band 
image of Gil. 2— 0.3 overlaid on the Chandra X-ray contours l IKoo et aI.II200'iT) . Boxes represent the locations of the three clumps observed 
in this study. 
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Fig. 2. — Integrated J- (a) and //-band (b) spectra of Clump 1, 
Clump 2, and Clump 3 over the entire region of each clump. The 
intensities are normalized by the peak [Fe II] 1.644 fim intensity of 
Clump 1. Constant offsets of 0.2 and 0.4 are subtracted from the 
spectra of Clump 2 and Clump 3, respectively, in order to show 
their spectra effectively. The spectra of Clump 3, which are faint, 
are multiplied by 5 for better display. In the bottom panel (c), the 
solid line is the observed [Fe II] 1.644 /im spectrum of SE-BHVF 
(see Figure [T] for the position) , whereas the dotted lines are two 
fitted gaussian components of -20 and -440 km s~^. 



Fig. 3. — Distribution of the [Fe II] 1.644 /^m line intensity 
(a), extinction (b), and electron density (c) estimated in Clump 1. 
Overlapped contours indicate [Fe II] 1.644 fim line intensities at 30, 
50, 70, and 90 % levels of the peak value. Vertical dotted line in (a) 
shows the positions along which we obtain one-dimensional radial 
profiles of the [Fe II] 1.644 fim line intensity and the [Fe II] 1.534 
to 1.644 fim line ratio (Figure|4]l. 
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Fig. 4. — One-dimensional radial profiles of the [Fe II] 1.644 fim 
intensity (filled circles) and the [Fe II] 1.534 to 1.644 ^m line ratio 
(empty circles) along the cut crossing the peak position of Clump 
1 (see Figure 3). The abscissa is the distance from the central 
pulsar of Gil. 2— 0.3 normalized by the distance to the shock front. 
The [Fe II] 1.644 fim line intensity is normalized by its maximum 
intensity. Also shown are model predicted profiles (see text for the 
details of the model): black-dotted line for [Fe II] 1.644 fim profile 
before beam dilution; red-solid line for the same [Fe II] 1.644 fim 
profile but after beam dilution; red-dashed line for beam-diluted 
[Fe II] 1.534 to 1.644 ^m line ratio profile. The green-solid line 
represents a beam-diluted [Fe II] 1.644 ^m profile from the model 
calculation of multiple filaments of finite sizes. 
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TABLE 1 

Measured line intensities of Clump 1, Clump 2, and Clump 

3 



Line 


Clump 1 


Clump 2 


Clump 3 


Transition 


A (/um) Observed Dereddened 


Observed Dereddened 


Observed Dereddened 



He I ^Si P 
[Fe II] ^ 0*^7/2 

[Fe II] ^ a^Dg/a 

[Fe II] o^Bt/.^, 
[Fe II] 061)3/2 ^ a^Dg/a 
[Fe II] a^F^i^ ^ a^Dg/g 
[Fell] a^Fy/a a* D3/2 
[Fc II] a-*F9/2 
[Fc II] a*Fr,/.2 
[Fc II] a-^FT/2 a^Dr,/2 
[Fe II] a*Fs/2 ^ a4D3/2 
[Fell] a'^F3/2^a*Di/2 
[Fe II] a^Fg/a ^ fl-^Dg/a 



1, 


.083 


3.5 ± 0.8 


50.0 ± 11.2 










1, 


.257 


30.3 ± 1.9 


136 ± 8.6 


25.9 ± 2.1 


136 ± 11.2 


22.5 ± 1.9 


136 ± 11.4 


1, 


.295 


3.4 ± 0.8 


12.7 ± 3.1 


3.5 ± 1.2 


14.7 ± 5.2 






1, 


.321 


15.2 ± 1.4 


49.4 ± 4.4 


11.2 ± 2.0 


41.0 ± 7.2 


9.8 ± 2.0 


40.0 ± 8.4 


1, 


.328 


5.7 ± 1.4 


17.8 ± 4.3 










1, 


.534 


13.7 ± 0.9 


19.2 ± 1.2 


12.8 ± 2.4 


18.5 ± 3.4 


8.5 ± 0.9 


12.7 ± 1.4 


1, 


.600 


6.1 ± 0.9 


7.0 ± 1.0 


9.0 ± 1.7 


10.4 ± 2.0 


8.5 ± 1.9 


9.9 ± 2.2 


1, 


.644 


100 ± 5.1 


100 ± 5.1 


100 ± 9.7 


100 ± 9.7 


100 ± 5.5 


100 ± 5.5 


1, 


.664 


4.4 ± 1.1 


4.1 ± 1.0 


4.7 ± 1.3 


4.4 ± 1.2 






1, 


.677 


8.9 ± 0.6 


8.2 ± 0.6 


11.6 ± 1.8 


10.5 ± 1.6 


9.3 ± 0.9 


8.3 ± 0.8 


1, 


.712 


3.4 ± 0.4 


2.8 ± 0.4 










1, 


.745 


3.8 ± 0.6 


2.9 ± 0.5 










1, 


.798 


6.0 ± 1.2 


4.1 ± 0.8 


9.7 ± 2.9 


6.4 ± 1.9 







Note. — Normalized to make the peak intensity of [Fe II] 1.644 
line of each clump 100. 



